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ABSTRACT 

The difference in the adiabatic gradient 71 between inverted solar data 
and solar models is analyzed. To obtain deeper insight into the issues of 
plasma physics, the so-called "intrinsic" difference in 71 is extracted, that is, 
the difference due to the change in the equation of state alone. Our method 
uses reference models based on two equations of state currently used in solar 
modeling, the Mihalas-Hummer-Dappen (MHD) equation of state, and the 
OPAL equation of state (developed at Livermore). Solar oscillation frequencies 
from the SOI/MDI instrument on board the SOHO spacecraft during its first 144 
days in operation are used. Our results confirm the existence of a subtle effect 
of the excited states in hydrogen that was previously studied only theoretically 
( [Nayfonov fc Dappenj 1998). The effect stems from internal partition function of 
hydrogen, as used in the MHD equation of state. Although it is a pure-hydrogen 
effect, it takes place in somewhat deeper layers of the Sun, where more than 
90% of hydrogen is ionized, and where the second ionization zone of helium 
is located. Therefore, the effect will have to be taken into account in reliable 
helioseismic determinations of the astrophysically relevant helium- abundance of 
the solar convection zone. 

Subject headings: helioseismology, equation of state, partition function, excited 
states 
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1. Introduction 

Helioseismology has proved to be a an extremely powerful tool to study the solar 
interior. Helioseismic inversion techniques have been used to probe details of the solar 
structure (cf. |Gough et al. 1996| | Basu et al. 1996j ; [Kosovichev et al. 1997] ), abundances 



(cf. Pasu fc Antia| ; [Kosovichev 19961 ; [Antia fc Chitre|) , and rotation rate (cf. [Thompson 



|et al. 19961 ; |Schou et al. 19981) . High-quality helioseismological data have allowed us for 



some time to go beyond studying solar structure and probe the properties of solar matter 
(see, e.g. ,fJhristensen-Dalsgaard et al. 1996|). Here, we try to study the equation of state of 



solar matter using inversions of solar oscillation frequencies. Our inversions are differential, 
that is, they determine the solar values in terms of deviations from given reference models. 
For the present analysis two sets of reference models have been used. The are based on 
the two equations of state normally used to construct solar models (see section 2.2), the 
so-called MHD equation of state flHummer fc Mihalas| 1988; [Mihalas et al. 1988| ; Pappcn 



et al. 198~8| ; |Dappen et al. 19871) , an d the OPAL equation of state ( [Rogers 1986| ; [Iglesias <S 



Rogers 1995 ; Rogers et al. 1996 and references therein), respectively. 



The observed solar oscillation modes are standing acoustic waves; hence the quantity 
most obviously probed is sound speed. Since the oscillations are largely adiabatic (except 
very near the surface), the frequencies are determined predominantly by the local adiabatic 
sound speed c, which is a thermodynamic quantity. In addition, the frequencies depend on 
the distribution of the density p in the Sun. Equivalently, adiabatic oscillation frequencies 
can also be determined from the distribution of density and the adiabatic gradient 71 in the 
Sun (see section 3.2 for more details). Since 71 is a dimensionless quantity, it is especially 
suited to study qualitative effects in the equation of state. In addition, for a high-precision 
diagnosis of the equation of state it is very fortunate that the Sun has a substantial 
convection zone, where equation of state effects show up independent of our uncertainty 
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in the opacity. The reason is that in the bulk of the convection zone, the stratification is 
essentially adiabatic and thus determined by thermodynamics ( |Christensen-Dalsgaard <S 
Dappen 199^ ). 



Helioseismic equation of state studies have shown the necessity to include the leading 
Coulomb correction (expressed by the Debye-Hiickel approximation) in order to meet the 
helioseismological constraints. This has been recognized in the early 80s and models with 



more sophisticated equations of state were tried in helioseismic studies (e.g. [Berthomieu et 



al. 198Q ; |Lubow et al. 1980| ; |Ulrich 1982| ; [Ulrich fc Rhodes 1983| ; [Shibahashi et al. 1983 



1984[ [Noels et al. 19841). With better data towards the end of the 80s, a clearer picture 



emerged ( |CJhristensen-Dalsgaard et al. 1988; Basu fc Antia 1995; Uhristensen-Dalsgaard fc 



Dappen 1992| ; |Christensen-Dalsgaard et al. 1996]) . However, even those equations of state 



that do include the leading Coulomb correction differ from the solar data by more than the 
observational uncertainty. They also differ discernibly among themselves, thus revealing the 
importance of the treatment of the nonideal effects beyond the leading Coulomb term. The 
aforementioned studies were based on the forward approach, i.e. , the comparison between 
observed and theoretically predicted solar oscillation frequencies. Equation-of-state studies 
based on inversions have also been attempted earlier (cf. pziembowski, Pamyatnykh fc 



Sienkiewicz 1992| ; |Elliot 19T?B| ; |Elliot fc Kosovichev| 1998), but they have been indirect. Those 



studies involved looking at either the sound-speed differences or differences in the adiabatic 
index 71, between the solar models and the Sun, both of which have contributions from 
differences in structure and composition in addition to contributions form the difference in 
the equation of state. 

In this paper we use the method proposed by Pasu fc Christensen-Dalsgaardl (1997) 
to study the difference in 71 between solar models that results from the differences in the 
equation of state alone, but not from the ensuing change in solar structure. To distinguish 
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it from the total 71 difference S^i, we call this contribution the "intrinsic" 71 difference 
(^7i)int- An immediate target of the new observational technique is the subtle effect of the 
excited states in hydrogen that that has so far only been studied theoretically ([Nayfonov 
fc Dappen] 1998). We begin with a presentation of the equation of state issues in order to 
discuss the hydrogen partition function effect. The solar oscillation frequencies used for this 
work were obtained from the SOI/MDI instrument on board the SOHO spacecraft during 
its first 144 days in operation ( jKhodes et al. 199"7| ). 



2. Equation of State Issues 

2.1. Ideal and Non-ideal Plasmas 

The simplest model is a mixture of nuclei and electrons, assumed fully ionized and 
obeying the classical perfect gas law. However, an "ideal-gas" equation of state can be 
more general. It may include deviations from the perfect gas law, namely ionization, 
radiation and degeneracy of electrons, as long as the underlying microphysics of these 
additional effects is still ideal, that is, as long as it does not contain interactions. The 
"particles", however, can be classical or quantum, material or photonic. In such an ideal 
framework, bound systems (molecules, atoms, ions) are allowed to have internal degrees 
of freedom (excited states, spin). All such ideal effects can be calculated as exactly as 
desired. Incidentally, in a recent study by [Elliot fc Kosovichev| (1998) it turned out that 
helioseismological accuracy demands inclusion of the relativistic effect in the electron gas, 
which has been neglected both in MHD and OPAL. Subsequently, both these equations of 
state are now being upgraded to include a relativistic treatment of electrons. Interestingly, 
the simple but astrophysically useful |Eggleton, Faulkner fc Flannery (1973)| equation of 



state (EFF) includes relativistic electrons. 
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One measure of nonideality in plasmas is the so-called coupling parameter T. In a 
plasma where particles have average distance (r) from each other, we can define T as the 
ratio of average potential binding energy over mean kinetic energy ksT (in the simplest 
case of hydrogen; generalizations to other elements are straightforward) 



Plasmas with T ^> 1 are strongly coupled, those with T 1 are weakly coupled. A famous 
example of a strongly coupled plasmas is the interior of white dwarfs, where the coupling 
can become so strong as to force crystallization. Weakly coupled are the interiors of stars 
with masses ranging from the slightly sub-solar ones to the largest. 

As one can suspect, T is the dimensionless coupling parameter according to which one 
can classify theories. Weakly-coupled plasmas lend to systematic perturbative ideas (e.g. 
in powers of T), strongly coupled plasma need more creative treatments. Improvements in 
the equation of state beyond the model of a mixture of ideal gases are difficult. This has 
both conceptual and technical reasons. As a fundamental conceptual reason we mention 
the fact that in a plasma environment already the idea of isolated atoms (and compound 
ions) has to be abandoned. A technical reason is the difficulty encountered when specific 
nonideal effect are modeled. The three principal nonideal effects are related to: (i) the 
internal partition functions of bound systems, (ii) pressure ionization, and (iii) collective 
interactions of the charged particles. The internal partition functions contain the difficult 
problem of excited states, where and how they are to be cut off. They are an important 
element in determining the ionization balances. Pressure ionization has to be provided by 
nonideal interaction terms, because ideal gases would unphysically recombine in the central 
regions of stars. 




k B T' 



(1) 
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2.2. Chemical and Physical Picture 

There are two basic approaches to realize nonideal equations of state: the so-called 
chemical and physical pictures. In the chemical picture one assumes that the notion of atoms 
and ions still makes sense, and ionization is treated like a chemical reaction. Modifications 
of atomic states by the surrounding plasma are expressed in a heuristic and intuitive way. 
A typical procedure is to find an expression for the probability that a given atomic state 
of an individual atom still exists. Examples of equations of state realized in the chemical 
picture are the MHD equation of state (see below) as well as the aforementioned simple 
EFF equation of state. 

The physical picture provides a systematic method to include nonideal effects. The 
approach starts out from the grand canonical ensemble of a system of the basic constituents 
(electrons and nuclei), interacting through the Coulomb potential. Configurations 
corresponding to bound combinations of electrons and nuclei, such as ions, atoms, 
and molecules, arise in this ensemble naturally as terms in cluster expansions. Any 
effects of the plasma environment on the internal states are obtained directly from the 
statistical-mechanical analysis, rather than by assertion as in the chemical picture. The 
only equation of state realized in the physical picture that has been applied to stellar 
interiors is the OPAL equation of state (see below). 

2.2.1. Chemical Picture 

Most realistic equations of state that have appeared in the last 30 years belong to 
the chemical picture and are based on the free-energy minimization method. This method 
uses approximate statistical mechanical models (for example the nonrelativistic electron 
gas, Debye-Hiickel theory for ionic species, hard-core atoms to simulate pressure ionization 



- 8- 



via configurational terms, quantum mechanical models of atoms in perturbed fields, etc.). 
From these models a macroscopic free energy is constructed as a function of temperature 
T, volume V, and the particle numbers Aq, . . . , N m of the m components of the plasma. 
This free energy is minimized subject to the stoichiometric constraint. The solution of this 
minimum problem then gives both the equilibrium concentrations and, if inserted in the 
free energy and its derivatives, the equation of state and the thermodynamic quantities. 

Obviously, this procedure automatically guarantees thermodynamic consistency. As an 
example, when the Coulomb pressure correction (to the ideal-gas contribution) is taken into 
account in the free energy (and not merely in the pressure), it affects both the pressure and 
the equilibrium concentration, i.e., the degrees of ionization. In contrast, the mere inclusion 
of the pressure correction would be inconsistent with other thermodynamic quantities. In 
the chemical picture, perturbed atoms must be introduced on a more-or-less ad-hoc basis to 
avoid the familiar divergence of internal partition functions (see e.g. |Ebeling et al.| 1976). In 
other words, the approximation of unperturbed atoms precludes the application of standard 
statistical mechanics, i.e. the attribution of a Boltzmann-factor to each atomic state. The 
conventional remedy of the chemical picture against this is a modification of the atomic 
states, e.g. by cutting off the highly excited states in function of density and temperature 
of the plasma. A currently popular equation of state realized in the chemical picture is 
the MHD equation of state ( [Hummer fc Mihalas| 1988; [Mihalas et al. 1988| ; |Dappen et al 



1988| ; Pappen et al. 1987|) , based on an occupation probability formalism. Specifically, the 



internal partition functions Zf 1 of species s adopted by MHD are weighted sums 



^ = £«;^«p(-^) . (2) 

Here,the subscript is labels states i of species s. Ei S are their energies, and the coefficients 
Wi S are the occupation probabilities that take into account charged and neutral surrounding 
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particles. In physical terms, Wi S gives the fraction of all particles of species s that 
can exist in state % with an electron bound to the atom or ion, and 1 — Wi S gives the 
fraction of those that are so heavily perturbed by nearby neighbors that their states are 
effectively destroyed. Perturbations by neutral particles are based on an excluded-volume 
treatment and perturbations by charges are calculated from a fit to a quantum-mechanical 
Stark-ionization theory (for details see |Hummer fc Mihalas] 1988). 



2.2.2. Physical Picture 

It is clear from the preceding subsection that the advantage of the chemical picture 
lies in the possibility to model complicated plasmas, and to obtain numerically smooth and 
consistent thermodynamic quantities. Nevertheless, the heuristic method of the separation 
of the atomic-physics problem from that of statistical mechanics is not satisfactory, and 
attempts have been made to avoid the concept of a perturbed atom in a plasma altogether. 
This has suggested an alternative description, the physical picture. In such an approach 
one expects that no assumptions about energy-level shifts or the convergence of internal 
partition functions have to be made. On the contrary, properties of energy levels and the 
partition functions should come out from the formalism. 

There is an impressive body of literature on the physical picture. Important sources 
of information with many references are the books by |Ebeling et al.| (1976), [Kraeft et al 



(1986), and [Ebeling et al.| (1991). However, the majority of work on the physical picture 



was not dedicated to the problem of obtaining a high-precision equation of state for stellar 
interiors. Such an attempt was made for the first time by the OPAL group at Livermore 
( [Rogers 198ft M 

esias fc Rogers 1995| ; [Rogers et al.| 1996 and references therein). 



To explain the advantages of this approach for partially ionized plasmas, it is 
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instructive to discuss the activity expansion for gaseous hydrogen. The interactions in this 
case are all short ranged and pressure is determined from a self-consistent solution of the 



equations flRogers 1981 



P 



knT 



z + z% + z% + 



p 



fdp\ 



(3) 
(4) 



k B T\dz) ' 

where z = A~ 3 exp(/i//c B T) is the activity, A = h/y / 27im e kBT is the thermal (de Broglie) 
wavelength of electrons, /i is the chemical potential and T is the temperature. The b n are 
cluster coefficients such that b<i includes all two particle states, 63 includes all three particle 
states, etc.. 

In contrast to the chemical picture, which is plagued by divergent partition functions, 
the physical picture has the power to avoid them altogether. An important example of such 
a fictitious divergence is that associated with the atomic partition function. This divergence 
is fictitious in the sense that the bound-state part of 62 is divergent but the scattering state 
part, which is omitted in the Saha approach, has a compensating divergence. Consequently 
the total 62 does not contain a divergence of this type ( [Ebeling et al.| 1976; [Rogers 19771) . 
A major advantage of the physical picture is that it incorporates this compensation at the 
outset. A further advantage is that no assumptions about energy- level shifts have to be 
made (see the previous subsection); it follows from the formalism that there are none. 

As a result, the Boltzmann sum appearing in the atomic (ionic) free energy is replaced 
with the so-called Planck-Larkin partition function (PLPF), given by (e.g. |Ebeling et al.| 
1976; [Kraeft et al.| 1986; [Rogers 1986| and references therein) 



PLPF = ^2(21 + 1) 



ni 



E, 



;l • 



(5) 
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The PLPF is convergent without additional cut-off criteria as are required in the chemical 
picture. We stress, however, that despite its name the PLPF is not a partition function, 



but merely an auxiliary term in a virial coefficient (see, for example, Dappen et al. 1987). 



2.3. Coulomb Correction and the Debye-Hiickel Approximation 

Debye-Hiickel (DH) theory is based on the replacement of the long-range Coulomb 
potential by a screened potential (see e.g. Lbeling et al.| 1976). This interaction leads 



to a negative pressure correction to the ideal-gas value. Under solar conditions, the 
relative correction culminates at about 8% in the outer part of the convection zone and 
it has another local maximum of about 1% in the core. Originally, the DH formalism 
included an additional ingredient of a fixed size for positive ions, inside of which the 
electrostatic potential is assumed to be constant. In astrophysical applications, such as 
MHD, this so-called r-correction has been adopted in the form given by Harris et al. (1960) 
(see also Gabriel 1994 ; Baturin et al. 1996| ). The systematic OPAL equation of state, 



however, comes close to the unmodified DH results under solar condition. This makes any 
r-correction rather implausible, especially since in comparisons with observational data, 
OPAL seems to fare better than MHD ( |Christensen-Dalsgaard et al. 1996|) , at least below 
the helium ionization zones where perhaps other effects, related to excited states, might 
dominate (see section 4). 

We would like to emphasize that close attention to the DH theory is warranted, 
because it describes the main truly nonideal effect under solar conditions. It was suggested 
in a number of early papers (e.g. Bcrthomieu et al. 1980 ; Ulrich 1982 ; Ulrich fc Rhodes 



1983| ; [Bhibahashi et al. 1983; 1984| ; [Noels et al. 19841) that improvements in the equation 



of state can reduce discrepancies between theory and observations. Later, |Christensen-| 
Dalsgaard et al.| (1988) showed that the MHD equation of state reduced these discrepancies 
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significantly for a large range of oscillation modes. It turned out that the Coulomb term 
is the dominant nonideal correction in the hydrogen and helium ionization zones. This 
discovery led to an upgrade of the simple, but astrophysically useful [Eggleton, Faulkner fc 



[Flannery (19 73 )| equation of state (EFF) through the inclusion of the Coulomb interaction 



term (CEFF) (|Christensen-Dalsgaard 1991| ; |Christensen-Dalsgaard fc Dappen 1992|) 



2.4. Beyond the Debye-Hiickel Correction 

Due to the relatively high temperature inside the Sun, the potential energy of the 
Coulomb interaction is small compared to the kinetic energy of particles, which allows 
us to believe that the DH theory makes quantitative sense at least asymptotically. But 
to estimate the possible error we need a physically based expression for next terms in 
the corresponding expansion. The OPAL equation of state contains such higher terms. 



An alternative is the path-integral based Feynman-Kac (FK) formalism of Alastueyfe 



Ferez| (1992, 1996), |Alastuey et al.| (1994, 1995). Its application to solar models is in 



progress ( |Perez fc Dappen 1998|) . 

New insight beyond DH might come through the equivalent description of the DH 
correction as a self energy of electrons and nuclei. A first study has revealed that the 
thermodynamic consequence of screened bound-state energies and a shifted continuum is 
hypersensitive to the details of the screening and the method with which the thermodynamic 
quantities are evaluated. No conclusion has been reached thus far but a preliminary study 
has shown that static screening alone comes close to the observational data, leaving little 
room for any thermodynamic influence of dynamic screening ( Arndt et al. 199"8| ). Screening 



has also been considered in connection with nuclear reaction rates, more specifically those 



responsible for the solar neutrino flux (e.g. Briiggen fc Gough 1997| , and references therein). 
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Very recently, [Nayfonov fc Dappen| (1998) examined the signature of the internal 
partition function in the equation of state. That study has revealed interesting features 
about excited states and their treatment in the equation of state. The MHD equation 
of state with its specific, density-dependent occupation probabilities (see section 2.2.1) is 
causing a characteristic "wiggle" in the thermodynamic quantities, most prominently in 
X p = (dlnp/dlnp)T, but equally present in the other thermodynamic quantities. 

EDITOR: PLACE FIGURE [I] HERE. 



Figures la and lb show the presence of excited states in ji for the case of a pure 
hydrogen plasma (panel b showing 71 relative to a simplified MHD equation of state 
containing only the ground-state contribution of hydrogen atoms). Temperatures and 
densities are taken from a solar model. Density is implied but not shown in the figure 
(for more details and different thermodynamic variables and chemical compositions, 



see [Nayfonov fc Dappenj 1998). Five cases were considered: (i) MHD [standard MHD 
occupation probabilities of [Hummer fc Mihalasj (1988)], (ii) MHDgs [standard MHD 
internal partition function of hydrogen but truncated to the ground state (GS) term], (iii) 
OPAL: OPAL tables [version of November 1996 of |Rogers etaT\ (1996)], (iv) MHD PL 
[MHD internal partition function of hydrogen, but replaced by the Planck-Larkin partition 
function, Eq. (|^), (v) MHD PL GS [MHD PL truncated to the ground state term]. The effect 
of the inclusion of the excited states in the internal partition function is manifest in the 
differences between MHD and MHDgs , and between MHDpl and MHDpl,gs; respectively. 
The effect of the different occupation probability of the ground and excited states shows up 
in the difference between MHD and MHDpl, and between MHDgs an d MHDpl gs- It was 
found that the presence of excited states is crucial. Also, the wiggle, which is a genuine 
neutral- hydrogen effect, is present despite the fact that most of hydrogen is already ionized. 
The qualitative picture does not change when helium is added ( [Nayfonov fc Dappenj 1998). 
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3. Inverting the Data 
3.1. Introduction 

We have performed differential inversions with respect to given reference models. We 
examine the influence of the equation of state by choosing two sets of models, one with the 
MHD and the other with the OPAL equation of state. Because of its differential nature, our 
inversion procedure is more reliable for reference models close to the real solar structure. 
Since diffusion has now become part of the standard solar model (|Christensen-Dalsgaard| 



et al. 1996|) , it is included in our reference models M1-M8 (see Table 1). For a study 
with artificial data, we have also considered a reference model without diffusion (M9). We 
use it to demonstrate that our inversion procedure applied to the theoretical frequencies 
of model M9 is able to reconstruct the intrinsic 71-difference with respect to reference 
model (M5). We found that this is indeed the case even for the relatively bad model M9, 
which is most likely farther away from the true solar structure than each of the reference 
models M1-M8 used for the real data. In addition we have tested the robustness of the 
inferred results despite uncertainties in solar model inputs by using a number of different 
solar models. Comparing the results of inversions that are made with respect to different 
reference models with the directly evaluated difference of the models themselves is a 
powerful tool to assess the quality of the inversions. 



3.2. Method 



An inversion for solar structure (e.g., Pziembowski, Pamyatnykh fc Sienkiewicz| 
19901 ; Pappen et al. 1991| ; |Antia fc Basu 1994| ; Pziembowski et al. 1994]) generally 



proceeds through a linearization of the equations of stellar oscillations around a known 
reference model. The differences between the structure of the Sun and the reference 
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model are then related to the differences in the frequencies of the Sun and the model by 
kernels. Non-adiabatic effects and other errors in modeling the surface layers give rise to 
frequency shifts ( Pox fc Kidman 1984j ; [Balmforth 1992| ; puenther 1994j ; [Rosenthal et al 



|1995| ; [Rosenthal et al. 1998|) which are not accounted for by the variational principle. In 



the absence of any reliable formulation, these effects have been taken into account in an ad 
hoc manner by including an arbitrary function of frequency in the variational formulation 
( pziembowski, Pamyatnykh fc Sienkiewicz 199U| ). Thus the fractional change in frequency 
of a mode can be expressed in terms of fractional changes in the structure of the model and 
a surface term. 

When the oscillation equation is linearized — under the assumption of hydrostatic 
equilibrium — the fractional change in the frequency can be related to the fractional 
changes in two of the functions that define the structure of the models. Thus, 



(cf. Pziembowski, Pamyatnykh fc Sienkiewicz 1990|) . Here Sui is the difference in the 



frequency Ui of the ith mode between the solar data and a reference model. The functions 
c and p are the sound speed and density. The kernels K\ „ and K\ n2 are known functions 
of the reference model which relate the changes in frequency to the changes in c 2 and p 
respectively; and Ei is the inertia of the mode, normalized by the photospheric amplitude 
of the displacement. The term -F sur f results from the near-surface errors. 

The conversion of the kernels for (c 2 , p) to those of the set [(71)^, u, Y], where, 
u = Pj p, P being the pressure and Y the helium abundance is discussed in [Basu fc 
CJhristcnsen-Dalsgaard| (1997). After the conversion, Eq. (^j) can be written as 
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int 



We use the technique of Optimally Localized Averages (OLA) to carry out the 
inversion. The purpose of the OLA technique QBackus fc Gilbert 1968t [Kosovichev 



1995|) is to construct appropriate linear combinations, J2i Ci(r )5uji/u!i, of the data with 
coefficients Cj(r ) so that the result represents a localized average of the quantity (£71/71) int 
at the radius r . This is possible since the (7i)mt averaging kernel J2iCi(r )K^,^ is a 
unimodular function localized at r = r , and the sums J2i c i( r o)Kl l Su/u, J2i CiirajKySY and 
T,iCiF sur{ (ui)/Ei are small. 

This objective is achieved by minimizing 



\ i / i,j 



with the constraint that the averaging kernel be unimodular, i.e., 



$>(r ) / 'K* Mi Jr)dr = 1 . (9) 

i 

Here, Eij is the covariance matrix of errors in the data. The parameters (3\ and 
control the contributions of 5u/u and SY, respectively, and fi is a trade-off parameter which 
controls the effect of data noise. The function w(r) is a suitably chosen, increasing function 
of radius, which ensures that the contributions from the second and third terms from the 
surface layers are suppressed properly. To reduce the influence of near-surface uncertainties 
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we apply the additional constraints that 

Ci^Er^fa) = , A = 0, . . . , A , (10) 

i 

where the $a are B-Splines with a suitably scaled argument (cf. Pappen et al. 1991] ). 

3.3. Reference Models 

We have constructed a set of calibrated solar models, all of which have either the 
MHD or the OPAL EOS. All models use OPAL opacities ([Iglesias fc Rogers 1996[) which 



are supplemented by the low temperature opacities of |Kuruc^ (1991). Since one of the 



most uncertain aspects of solar modeling is the formulation of the convective flux, we use 
two different convective flux formalisms — the standard mixing length theory (MLT) and 
that of |Canuto fc Mazitelli| (1991) (henceforth called CM). The two formalisms give fairly 



different stratifications in the outer regions of the Sun where (7i)i n t differences are also 
likely to be the largest. All models have Z/X = 0.0245 (cf. £3 revesse fc Noels| 1993). 

EDITOR: PLACE TABLE | HERE. 

A major source of uncertainty in the inversion results is the solar radius. Till recently, 
the standard value of 695.99 Mm was used unquestioningly. However, new f-mode data 
suggests that the radius actually smaller (cf. |5chou et al. 1997] [Antia 19981) . A difference 



in the radius between the solar models used as a reference model for inversions and the Sun 
can lead to a large error in the inversion results (cf. Pasu 1998a]) . Thus we have constructed 



models with the standard radius as well as the f-mode radius of 695.78 Mm (cf. |Antia 19T?5| ). 
Recent observations to determine the solar radius (cf. [Brown fc Christensen-Dalsgaard 1998| ) 
suggest that the solar radius is even lower, 695.508 Mm. Some of the relevant properties of 
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the models are listed in Table [l|. All models, except M9, assume the gravitational settling 
of helium and heavy elements. The composition profiles for models M1-M8 are the ones 
obtained from helioseismological inversions by [Antia fc Chitrej (1998). For model M9, the 



profile from the no-diffusion model of [Bahcall fc PinsonneauTt| (1992) was used. 



Note that the surface helium abundances and convection zone depths of all models 
except M9 are consistent with the helioseismically determined helium abundance and depth 
of the convection zone of the solar envelope (cf. | Basu fc Antia| 1997; [Basu 1998a]) . 



3.4. Test results 

To test the inversion procedure, we have inverted models Ml and M9 using model M5 
as the reference model. The purpose of inverting model Ml is to show that we can invert 
for the intrinsic 71 differences between two models (Fig. 2a). We see that we are indeed 
successful in inverting for the intrinsic 71 difference between MHD and OPAL models. The 
structure and helium abundance of the models Ml and M5 are very similar, hence the the 
71 difference between models with different equations of state is dominated by differences 
in the EOS rather than differences in structure. 



EDITOR: PLACE FIGURE @ HERE. 



In a parallel test, model M9 is inverted with respect to reference model M5 in order to 
show that we are indeed inverting just for the equation of state part of the 71-difference 
(^7i)int- We deliberately choose for this test the relatively bad solar model M9 without 
diffusion, to show insensitivity to the quality of the reference model. Importantly for 
this test, model M9 has the same equation of state as the reference model M5. Using 
theoretically computed frequencies of model M9, our inversion procedure gives the intrinsic 
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difference of 71 between models M5 and M9. Since the equation of state is the same in the 
two models, we should obtain essentially zero for the intrinsic difference. This is indeed the 
case (see Fig 2b). 

The small deviations from zero can be explained mainly as contributions in the solution 
because of the large difference in the helium abundance [i.e. the third term in Eq. (0)], 
though there can also be a contribution from the neglected differences in the heavy-element 
abundances. Increasing the parameter @ 2 helps in getting better inversion results, but it 
degrades the quality of the inversion. However, this is not a concern for inversion of the 
solar data with models M1-M8, since the difference in the helium abundance between these 
models and the Sun is an order of magnitude less than the helium-abundance difference 
between models M9 and M5. For comparison, Fig. 2b also shows total 71 difference. Note 
that the total 71 inversion captures the difference in total 71 to within the limits of the 
resolution. 



EDITOR: PLACE FIGURE E HERE. 



The mode set and errors used for the test were the same as that of the SOI/MDI set 
mentioned above ( [Rhodes et al. 1997|) . We are quite successful in constructing localized 
kernels till fairly close to the solar surface. A sample of the averaging kernels is shown 
in Fig. 3. The use of the observational mode-set and errors ensures that the inversion 
properties (e.g. resolution of the averaging kernels, propagated errors etc.) will remain the 
same when solar data are inverted. 
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4. Results and Conclusion 

The result obtained by inverting solar data using models M1-M8 are shown in Fig. 4. 
We show inversions with respect to a deliberately large number of reference models, having 
different parameters. By doing so we demonstrate that in the chosen technique only the 
equation of state matters. 

Note that the results of all the MHD models are very similar despite differences in 
radius and stratification. However these are very different from the OPAL models which 
form a distinct group. The differences between the two equations of state and the Sun are 
fairly large in the outer layers of the Sun 



EDITOR: PLACE FIGURE g HERE. 



In contrast to earlier results ( Pasu fc Christensen-Dalsgaard| 1997), which could only 
resolve the somewhat deeper layers of the Sun, the present focus on the 20% uppermost 
layers confirms the aforementioned effect of excited states as included in the MHD equation 
of state (section 2.4). The results also confirm that the OPAL equation of state is better 
for r < 0.97i? Q ( |Christensen-Dalsgaard et al. 1996| ), but the situation is reversed in the top 
2-3 per cent of the radius. 

Since the difference in ji between MHD and OPAL is the wiggle of Fig. 1, the observed 
preference of the MHD model in the upper region could indicate the correctness of an 
MHD-like treatment of the exited states. Again, we emphasize that Fig. 4 confirms earlier 
results that below the wiggle region, OPAL fares better than MHD ( |Christensen-Dalsgaard 



et al. 1996|) . There, the internal partition functions of OPAL a la Planck-Larkin, and the 



absence of a r-correction (section 2.3) seem to be the better choice. 

The reversal of fortune in favor of MHD in the upper part of the Sun (above 0.97 solar 
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radii) could be due to the different implementations of many-body interactions in the two 
formalisms. Since density decreases in the upper part, OPAL by its nature of a systematic 
expansion, inevitably becomes more accurate; but MHD might, by its heuristic approach 
(and by luck!), have incorporated even finer, higher-order effects. 

The helioseismic helium-abundance determination procedure is quite sensitive to the 
details of internal partition functions. Since the wiggle is located in the second ionization 
zone helium, although it is a pure hydrogen effect (see section 2.4), it has a very important 
bearing on the astrophysically relevant helioseismic helium-abundance Y determination in 
the solar convection zone. In principle, the asymptotic analysis by |Basu fc Antia| (1995), 
which leads to Y — 0.249 for OPAL and Y = 0.246 for MHD, could be taken as indication 
that the equation of state does not matter in the result for Y. However, other studies came 
to a different conclusion. In independent inversions, [Kosovichev| (1997) obtained = 0.23 for 
MHD and Y = 0.25 for OPAL, |Bichard et al| (1998) Y = 0.242 for MHD and Y = 0.248 
for OPAL, and gagg (1998b) Y = 0.248 for OPAL and Y = 0.251 for MHD. In an entropy 
calibration by |Baturin fc Ayukov| (1997) the result was Y = 0.23 for OPAL and Y = 0.25 
for MHD. In the latter two, the sign of the change is reversed compared to the three former. 

The confirmation of the excited-states effect in solar data makes helium-abundance 
determinations based on the MHD equation of state more plausible than the ones based 
on OPAL, but the various current techniques do not yet allow a consensus. We note in 
passing that the helium- abundance values given in Table 1 should not be confused with 
the aforementioned inversion results. Table 1 lists the result of calibrations in which the 
helium parameter compensates for other effects as well, and the relative stability of the 
helium parameter in M1-M8 is no indication for its absolute correctness. Inversions can 
in principle isolate the helium abundance from other effects, but an absolutely accurate 
equation of state is required. As we mentioned, any quantitative conclusions about the 
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value of Y would be premature. However, since the MHD equation of state appears to 
describe reality in the helium-ionization zone better, it might be the preferable equation of 
state in helium-abundance determinations. 

Let us add a word of caution, though. It could appear tempting to produce a 
"combined" solar equation of state, with MHD for the top part and OPAL for the lower 
part. However, such a hybrid solution is fraught with danger. For instance, it is known that 
patching together equations of state can introduce spurious effects flDappen et al. 1993|) . 
It seems that the right way is to improve MHD and OPAL in parallel and independently, 
guided by the progress of helioseismology. 
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FIGURE CAPTIONS 

Fig. 1. — Panel (a): absolute values of 71 for solar temperatures and densities of a hydrogen- 
only plasma. Linestyles: MHD - asterisks, MHD G s - dashed lines, MHD PL - dotted-dashed 
lines, MHDpl gs ~~ dotted lines, and OPAL - solid lines. Panel (b): relative differences with 
respect to MHD G s, in the sense (71 — 7i[MHD G s])/7i[MHD G s]), using the same line styles 
as in (a). The horizontal solid zero line, representing MHD GS , is also shown. See text for 
the definitions of the different MHD versions. 

Fig. 2. — Panel (a): Result of inverting model Ml using model M5 as the reference model. 
The continuous line is the exact difference in the intrinsic 71 between the two models. The 
points are the results obtained by inverting the frequency differences between the two models. 
The vertical error bars represent the errors in the solution because of observational errors in 
the frequencies. The horizontal error bars are a measure of the resolution of the inversion 
and is the distance between the quartile points of the averaging kernels. Panel (b): Result 
of inverting model M9 with model M5 to obtain the intrinsic and total 71 difference. The 
continuous line is the exact difference in total 71 between the two models. The models were 
constructed with the same EOS, hence the intrinsic 71 difference is zero. 

Fig. 3. — Averaging kernels obtained at selected radii. 

Fig. 4. — Relative difference between 7! obtained from an inversion of helioseismological 
data and 71 of the solar models M1-M8 listed in Table I, in the sense "Sun - model". Only 
the "intrinsic" difference in 71 is shown, that is, the part of the difference due to the equation 
of state (see text). The filled points are results obtained with MHD models, the empty ones 
are results with OPAL models. Lines have been drawn through results of models Ml and 
M5 to guide the eye. Error bars have been drawn only on two sets of results for the sake of 
clarity. 
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Table 1: Properties of the solar models used in Figs. 2a, 2b, and 4 (see text). 



Model EOS Radius Convective Y s r cz / R Q 







Mm 


Flux 






Ml 


MHD 


695.78 


CM 


0.2472 


0.7145 


M2 


MHD 


695.99 


CM 


0.2472 


0.7146 


M3 


MHD 


695.51 


CM 


0.2472 


0.7145 


M4 


MHD 


695.78 


MLT 


0.2472 


0.7146 


M5 


OPAL 


695.78 


CM 


0.2465 


0.7134 


M6 


OPAL 


695.99 


CM 


0.2465 


0.7135 


M7 


OPAL 


695.51 


CM 


0.2466 


0.7133 


M8 


OPAL 


695.78 


MLT 


0.2465 


0.7135 


M9 


OPAL 


695.78 


CM 


0.2646 


0.7268 



